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I. INTRODUCTION 

New functionalities can be realized by integrating fer- 
romagnetic elements into electronic circuits and devices. 
The interplay between magnetism and electric currents 
in these structures is utilized by the giant magnetoresis- 
tance (GMR), the operating principle of the read heads in 
modern magnetic hard disk drives. Considerable progress 
has been made in improving magnetic random access 
memories.^ Efforts to further miniaturize and improve 
the performance of magnetoelectronic devices are ongo- 
ing in academic and corporate laboratories. Low power 
consumption and noise levels are essential. In spite of the 
technological relevance, a comprehensive understanding 
of coupled current and magnetization noise and the re- 
lated energy dissipation in nanoscale magnetoelectronic 
circuits is lacking. 

From the early studies of Johnson^ and Nyquistji we 
know that the equilibrium voltage noise power in con- 
ductors is proportional to the electric resistance. This 
relation between the equilibrium noise and the out-of- 
equilibrium energy dissipation is a standard example of 
the fluctuation-dissipation theorem (FDT).*'^ In recent 
years, important advances have been made in the un- 
derstanding of electronic equilibrium (thermal) and non- 
equilibrium (shot) noise in mesoscopic conductors.^ 

The electron spin plays an important role in elec- 
trical noise phenomena in magnetic multilayers. In 
early theoretical studies^i'^'^ i-'^^'^^i^^ of charge and spin- 
polarized current noise in such systems, magnetizations 
were assumed to be static. However, the magnetiza- 
tion itself fluctuates as well. Thermal fluctuations of 
the magnetization vector in isolated single-domain fer- 
romagnets have been analyzed by Brown, ^"^ who intro- 
duced a stochastic Langevin field acting on the magne- 
tization to account for thermal agitation. His proof that 
this field's (white-noise) correlator is proportional to the 
magnetization damping (see below) is another manifes- 
tation of the FDTfiiii^ The stochastic field can be in- 



troduced into the spatiotemporal equation of motion for 
the magnetization (Landau-Lifshitz-Gilbert equation), 
affecting, e.g., current-driven magnetization dynamics 
and reversaliiSiilii^ii^ 

A moving magnetization vector in ferromagnets un- 
dergoes viscous damping that relaxes the magnetization 
toward the lowest (free-)energy configuration. This pro- 
cess is in practice well described by a phenomenological 
damping constant, introduced by Gilbert. Despite 
some progresS j^^i^^'^'^i^^'^^i^^ a rigorous quantitative un- 
derstanding of the magnetic damping in transition-metal 
ferromagnets has not yet been achieved. The theory of 
the enhanced Gilbert damping in ferromagnets in good 
electrical contact with a conducting environment is in a 
better shape. The loss of angular momentum due to spin 
current pumping into the environment agrees with the 
Gilbert phenomenology, ■^^'^^ and experiment and theory 
addressing the additional damping agree well with each 
other. 

The electronic and magnetic fluctuations in magne- 
toelectronic structures are intimately coupled to each 
otheri^Si^ For example, the magnetization noise in fer- 
romagnetic films in good electric contact with normal 
metals has been predicted to increase due to spin cur- 
rent fluctuations: The spin current components polar- 
ized perpendicularly to the magnetization are absorbed 
at the interface, leading to a fluctuating spin-transfer 
torqueHi^i^i^ that induces additional magnetization 
noise. This noise is related to the excess Gilbert damp- 
ing caused by the angular momentum loss due to spin 
pumping, in accordance with the FDT. 

Here we investigate the interplay between (zero fre- 
quency) current and magnetization noise in multilayers 
of alternating magnetic and non-magnetic films. We 
take advantage of the FDT to relate the equilibrium 
electric (current and voltage) and magnetic (magneti- 
zation and field) noise to the corresponding dissipation 
of energy. We start by reviewing the noise in a single 
monodomain ferromagnet sandwiched by normal met- 
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als, including technical details that were omitted in Ref. 
Isol . Both thermal equilibrium (Johnson-Nyquist) current 
noise and nonequilibrium shot noise are taken into ac- 
count. Next, we consider spin valves, i.e., two ferromag- 
netic films separated by a normal metal spacer."^® We 
consider both a symmetric structure in which both lay- 
ers fluctuate, as well as an asymmetric one, in which 
one layer is assumed fixed. Magnetoelectronic circuit 
theory^i^^i^ is used to calculate the charge and spin 
current fluctuations. The resulting enhanced magnetiza- 
tion noise and Gilbert damping in principle are tensors 
that depend on the magnetic configuration. 

Spin valves provide an opportunity to indirectly mea- 
sure magnetization noise via resistance fluctuations, 
which are manifested by voltage noise for a current-biased 
system or current noise for a voltage-biased system.— 
This offers an experimental test of our theory. We ob- 
tain analytical expressions for the magnetic contribution 
to the induced electric noise for different magnetic con- 
figurations. The noise is of potential importance for the 
performance of spin valve read heads 4i For symmetric 
structures in which both layers fluctuate, dynamic cross 
talk between the layers becomes important, causing a 
possibly large difference in noise level between the paral- 
lel and antiparallel magnetic configurations. Our results 
for these spin valves include previously presented findings 
as a limiting case.^^^ After the completion of this work,"'^ 
it was shown that spin- valves in equilibrium also exhibit 
colored voltage fluctuations caused by spin pumping of 
the moving magnetizations.'*'^ 

The paper is organized as follows. We begin by review- 
ing the fluctuation-dissipation theorem, applied to mag- 
netic systems. In section III, the noise properties of a sin- 
gle ferromagnetic thin film sandwiched by normal metals 
is worked out in detail, emphasizing the relation of the 
noise to the damping. In section IV, we consider current 
noise, magnetization noise, and magnetization damping 
in spin valves, and use the results to calculate the resis- 
tance noise induced by GMR. In section V we summarize 
our conclusions. 



II. FLUCTUATION-DISSIPATION THEOREM 

The fiuctuation-dissipation theorem (FDT) relates the 
spontaneous time-dependent changes of an observable of 
a given system in thermal equilibrium its linear response 
to an external perturbation that couples to that observ- 
able. For example, in an electric conductor the sponta- 
neous fluctuations in the electric current are proportional 
to the dissipative (real) part of the conductivity, i.e., the 
response function to an applied electric field. ^'^ Similarly, 
the equilibrium fluctuations of the magnetization vector 
in a ferromagnet are proportional to the dissipative part 
of the magnetic susceptibility, i.e., imaginary part of the 
response function to an applied magnetic fleld. In the 
following, we briefly recapitulate this FDT for magnetic 
systems. 



Sufficiently below the Curie temperature, changes in 
the modulus of the magnetization are energetically costly 
and may be disregarded. For sufficiently small magnetic 
structures spin waves freeze out of the problem. Hence, 
a small ferromagnetic particle or thin film is well de- 
scribed in terms of a single magnetization vector Mgrn, 
where Ms is the magnitude of the magnetization and m 
a unit vector ( "macrospin" model) . The time-dependent 
equilibrium fluctuations of the magnetization are charac- 
terized by the autocorrelation function {Smi(t)Smj{t')), 
where 5mi(t) — •mi{t) — {rriiit)) are transverse fluctua- 
tions. Here the brackets denote statistical averaging at 
equilibrium, and i and j denote Cartesian components 
perpendicular to the equilibrium/average magnetization 
direction. The classical FDT states that these fluctua- 
tions are related to the magnetic susceptibility: 



{5mi{t)5mj{t')) 



'2t:MsV 



(1) 



where T is the temperature, V the volume of the ferro- 
magnet, and Xiji^) the ij-component of the transverse 
magnetic susceptibility at frequency uj. The latter is 
the linear (causal) response function that describes the 
changes of the magnetization, Ami(i), caused by an ex- 
ternal driving fleld H('^'')(i): 



Am,(t) 



j 



dt'xtj (t 



(2) 



An alternative form of the FDT that turns out useful in 
the course of this paper can be derived by introducing 
a stochastic magnetic field h^"^ (t) with zero mean. This 
field effectively represents the coupling of the magnetiza- 
tion to the dissipative degrees of freedom, and is viewed 
as the cause of the thermal fiuctuations Sm{t). The mi- 
croscopic origin of h^"^ (t) does not concern us here, but it 
might, e.g., represent thermally excited phonons that de- 
form the crystal anisotropy fields. From Eq. ^ it follows 
that Smi{LL!) = J2j Xiji^)hf'\Ltj) in frequency domain. 
Inverting this relation, the correlator of the stochastic 
field has to obey the relation 



2'kMsV j 



(3) 



where Xij^ ('^) is the zj-component of the Fourier trans- 
formed inverse susceptibility. 



III. SINGLE FERROMAGNET 

The magnetization dynamics of an isolated single- 
domain ferromagnet is well described by the Landau- 
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Lifshitz- Gilbert (LLG) equatioi 
dm 

-J- = -70m X Hcff - 
at 



.20.44 



dm 



aom X (4) 



where 70 is the gyromagnetic ratio, Hoff the effective 
magnetic field, and ag the Gilbert damping constant. 
The effective field has contributions due to crystal and 
form anisotropies, as well as externally applied magnetic 
fields. By linearizing this LLG equation we can evaluate 
the magnetic susceptibility and the equilibrium magne- 
tization noise. The average equilibrium direction of the 
magnetization is aligned with Hcff to minimize the en- 
ergy: mo = Hcff /I Hcff I . A weak external driving field 
is included by substituting Hcff Hcff + H^'*'^^(/;). In 
the present model only the component of H'"^''^ trans- 
verse to the magnetization will solicit a response m{t) « 
mo -I- Am(i) of the magnetization. Here Am{t) is nor- 
mal to mp. To lowest order in Am{t), the LLG equation 
gives the inverse susceptibility tensor matrix 



X 



1 

70 



70 |Hcff| - iwao iuj 

-iuj 70 |Hcff I - iujaa 



(5) 



in the plane normal to mo. A dependence of the effective 
field Hcff on m does not affect the noise properties. 

The magnetization noise follows from substituting 
Eq. ([5]) into Eq. ([T]) . The correlator of the stochastic field 
is obtained from Eqs. ^ and (O and does not depend 
on the effective fieldJ^ 



{hf\t)h 



"0 



6,,d{t - 1'). (6) 



The relation between the equilibrium magnetization fluc- 
tuations and the dissipation in the form of the Gilbert 
damping is evident. 

Up to now we considered a ferromagnet isolated from 
the outside world. Its dynamics is altered by embed- 
ding into a conducting environment A ferromagnet 
with time-dependent magnetization "pumps" an angular- 
momentum (spin) current 



' 47r V dt 



Img 



n 



dm 

~dr 



(7) 



into an adjacent conductor. Here g^^ is the dimension- 
less transverse spin ("spin mixing" ) conductance that 
depends on the interface transparency between ferromag- 
net and proximate metah ^^'^^ i^ When the spin current 
is efficiently dissipated in the conductor, thus does not 
build up a spin accumulation close to the interface, the 
loss of angular momentum corresponds to an extra torque 
7lP""^P/(M,y) on the right hand side of the Eq. Q. 
This is equivalent to an increased Gilbert damping and 
a modified gyromagnetic ratio:^- 
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7o^Img 
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7 f , Tofi-Reg 
ao ^ a = — ao + 

70 V 4:TTMsV 



(8) 
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In the strong coupling limit (intermetallic interfaces), 
Img^^ ^ Reg^^ and we are allowed to disregard the dif- 
ference between 7 and 70. 

Another term which modifies the magnetization dy- 
namics is the so-called spin-transfer torque i^i^^i^i^^ 
It is also proportional to the spin-mixing conduc- 
tance introduced above^^i^^i^ and represented by adding 
— 7oIs,abs/ {MsV) to the right hand side of Eq. (ji]). Here 
Is,abs is the spin-polarized current transversely polarized 
to the magnetization, which is absorbed by the ferromag- 
net on an atomic length scale, thereby transferring its 
angular momentum to the magnetization. Spin pump- 
ing and spin-transfer torque are related by an Onsager 
reciprocity relation4^ 

Recently we have shown^*' that the magnetization noise 
in magnetoelectronic nanostructures can be considerably 
increased as compared to an isolated ferromagnet. At 
elevated temperatures, thermal fluctuations in the spin 
current exert a fluctuating torque on the magnetization, 
increasing the noise. For a ferromagnet sandwiched by 
normal metals, the enhancement of the noise is described 
by a stochastic field h''^^ (t) similar to the intrinsic field 
h'"'(<). Its correlation function reads^ 



(/if)(OC^(i')) 



2fcsT 
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where 



n 



jhReg 
in MsV 



S^j6{t-t'), (10) 
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is the enhancement of the Gilbert damping due to spin 
pumping (see Eq. Assuming that h'^°)(t) and 

h'^'^)(i) are statistically independent, the total magne- 
tization noise is thus given by h{t) = h(°)(t) -I- h'-*'^\t). 
We know that the total damping is determined by a = 
aQ + a', and from Eqs. ^ and PTI| we see that the total 
noise is related to the total damping, in agreement with 
the FDT. Hence, the thermal spin current noise is the 
stochastic process related to the enhanced dissipation of 
energy by spin pumping. By calculating the noise power 
we also know the damping and vice versa. In thin ferro- 
magnetic films, a' can be of the same order or even larger 
than Oq.^^ In the following subsections we will give a de- 
tailed derivation of Eq. (fTU|) . We also evaluate the shot 
noise contribution to the magnetization noise, which is 
important at low temperatures We note here that Eq. 
PH)) may be found also by direct application of Eq. ([3]) 
to the LLG equation with spin-pumping included. 



A. Scattering theory 

We study a thin ferromagnetic film connected to two 
normal reservoirs, as shown in Fig. [1] The reservoirs 
are perfect spin sinks and the ferromagnet is taken to 
be thicker than the magnetic coherence length Ac = 
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FIG. 1: A thin ferromagnetic (F) film is sandwiched by nor- 
mal metals (N). The current fluctuations in the system are 
evaluated in terms of transmission probabilities for the elec- 
tron states, with the aid of second quantized annihilation 
and creation operators. The operators shown in the figure 
are annihilation operators, with the o-operators annihilat- 
ing electrons moving towards the ferromagnet, and the b- 
operators annihilating electrons moving away from the fer- 
romagnet. Also shown are the refiection and transmission 
matrices r,r' ,t,t' (see Eq. (|14[) '). for simplicity without spin 
indices. 



7r/(fc| — fcj), where are spin-dependent Fermi mo- 

menta. For transition metals, Ac is of the order of mono- 
layers. The normal metals are characterized by Fermi- 
Dirac distribution functions /l and fn with chemical 
potentials and ^r, where L and R refer to the left 
and right sides at a common temperature T. We use the 
Landauer-Biittiker (LB) scattering theory^ to evaluate 
the spin current fluctuations, and the LLG equation to 
calculate the resulting magnetization noise. 

In the LB approach electron transport is expressed in 
terms of transmission probabilities between the electron 
states on different sides of a scattering region. Here we 
interpret the ferromagnetic film as a scatterer that limits 
the propagation of electrons between the normal reser- 
voirs. The scattering properties of the ferromagnet and 
the bias between the reservoirs determine the transport 
properties of the system. The transport channels in the 
leads are modelled as ideal electron wave guides in which 
the transverse and longitudinal motions are separable. 
The transport channels at a given energy E are then la- 
beled by the discrete mode index for the quantized trans- 
verse motion, by which the continuous wave vector for the 
longitudinal motion is fixed. The LB formalismSi^ gen- 
eralized to describe spin transport leads to the current 
operator 



If it) 



h 



X [a\p{E)aAc.{E') ~ b\^,{E)bAa{E% (12) 



at time t on side A [— L(left) or i?(right)] of the fer- 
romagnetic film. Here, a and (3 denote components in 
2x2 spin space. aAa{E) and hAa{E) are operators for 
all transport channels at energy E that annihilate elec- 
trons with spin a in lead A that move towards and away 
from the ferromagnet, respectively (see Fig. [T|). The a- 



operators are related to the ^-operators by the scattering 
properties of the ferromagnet: 

bA^{E)^Y. 

SABaf3{E)aBf3{E), (13) 

BI3 

where is the scattering matrix for incoming elec- 

trons with spin /3 in lead B (= L oi R) scattered to 
outgoing states in lead A with spin a. The summation is 
over B = L, R and over spin /3 =f , j. A similar relation 
holds for the creation operators. Current conservation 
implies that the scattering matrix is unitary. Suppress- 
ing spin indices for simplicity (see Fig. [1]) 



bL 
bR 



r t' 
t r' 



aL 
aR 



(14) 



where r = sll, r' = srr, t = srl and t' = slr- In 
the following we disregard spin-flip processes in the fer- 
romagnet. Choosing the spin quantization z-axis in the 
direction of the average magnetization, this implies that 

SABaP = SABa^ap- 

The outgoing charge and spin currents are given 
respectively by Ic,A{t) = Y^a^Tit) and = 

-(?i/2e)X]Q/3^"^-^A"(*)' where a = (ct^,, CTj^, a^) is 
the vector of Pauli matrices. The expectation val- 
ues for charge and spin currents are evaluated using 
the quantum statistical average {a\,-^^{E)aBnf3{E')) = 
SABSmnSaf3S{E — E')fA{E) of the product of one creation 
and one annihilation operator, where m and n label the 
transport channels. The creation and annihilation oper- 
ators obey the anticommutation relation 

WAmJE),aBnpiE')} = SABSmnS^pSiE ~ E'), (15) 

whereas the anticommutators of two creation or two an- 
nihilation operators vanish. Similar relations hold for the 
b operators. The average 

(«Afea {Ei)aBip iE2)al^^ {E:i)aDnS {E4)} 

- i^Akc {Ei)aBip {E2)){ah^^ {E^)aDn5 [E,)) 

= SADSBC^krArnSasSp-ySiEi - Ei)5{E2 - E3) 

X fAiE,)[l - fB{E2)], (16) 

where the subscripts A,B,C,D denote leads, k,l,m,n 
transport channels, and a,/3,7,(5 spin, is needed for the 
calculation of the current fluctuations. We also need the 
identity 



E 

CD 



T^'^{sACa^ADps''BDl3^BCa) = SabMa , 



(17) 



which follows from the unitarity of the scattering matrix. 
Here the trace is over the space of the transport channels, 
and Ma is the number of transverse channels in lead A, 
all at a given energy. 

The charge and spin current correlation functions read 



,AB 



{t-t')^{5h,A{mc.B{t')) 



(18) 
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and 

S^J,AB{t - t') = {6h,AtWs,,B{t')), (19) 

where 5Ic.A{i) = Ic.A{t) — {Ic.A{t)) denotes the deviation 
of the charge current from its average value in lead A 
at time t, and SIs^.Ait) is the deviation of the vector 
component i (i = x,y or z) of the spin current. We 
are interested mainly in the low- frequency noise, i.e., the 
time integrated value of the correlation functions: 

Sc,ab{^ = 0) - j d{t~ t')S,^AB{t - t'). (20) 

Two fundamentally different types of current noise have 
to be distinguished: Thermal (equilibrium) noise and 
(non-equilibrium) shot noise. In general, the total noise 
is not simply a linear combination of both types. Nev- 
ertheless, it is convenient to treat the two noise sources 
independently, by separately investigating the noise of an 
unbiased system at finite temperatures in Sec. IIII Bl and 
the shot noise under an applied bias at zero temperature 
in Sec. [mCl 

B. Thermal current noise 



current 

CiA = ^WgT-,5i) (23) 

differs from the charge current noise only by the squared 
conversion factor, (ft/2e)^, from charge to spin currents. 
The transverse (polarized perpendicular to the magneti- 
zation) spin-current components fluctuate as 

S^!aa = S^^^AA = ^ksTign + ,iT). (34) 

The "spin mixing" conductances g]j = Tr[l — ftl'^i)^] — 
{g[^)* and = Tr[l - r\{r[y] = {g\^y parametrize 
the absorbtivity of the ferromagnetic interfaces for 
transverse-polarized spin currents. We see that also the 
spin-current noise obeys the FDT, since the spin-current 
correlators are proportional to the conductances for the 
respective spin current components. 

The cross correlation S'^*'']^^ — ~S^^\j^ reflects con- 
servation of the longitudinal spin current in the ferro- 
magnet, since spin-flip scattering is disregarded. On the 
other hand, S^^^^j^ — Syy^^^j^ — 0, because the transverse 
spin current is absorbed at the interfaces to a ferromag- 
net thicker than the magnetic coherence length. 



At equilibrium ~ fn ~ f , and the average current 
vanishes. However, at finite temperatures, the occupa- 
tion numbers of the electron channels incident on the 
sample fluctuate in time and so does the current. Using 
Eqs. dni), dlSl), dUl), dni) and /(l-/) = kBT{-df/dE), 
we recover the well-known Johnson-Nyquist noise 

sl^SAiu; = 0) = ^kBT{g^+g^) (21) 

in the zero-frequency limit. Here g" — Tr(l — rl^Ta), 
where the trace indicates again a summation over trans- 
port channels, is the spin-dependent dimensionless con- 
ductance of the ferromagnet, to be evaluated at the Fermi 
energy. The superscript (th) emphasizes that the fluctu- 
ations are caused by thermal agitation. The result for 

^'cAsi^ ~ 0), where B A, differs from the above ex- 
pression only by a minus sign, since current direction is 
defined positive towards the ferromagnet on both sides, 
and charge current is conserved. The Johnson-Nyquist 
noise, Eq. pT|) . is a manifestation of the FDT, since it 
relates the equilibrium current noise to the dissipation of 
energy prarameterized by the conductance. 

The thermal spin current noise can be obtained in a 
similar way. At zero frequency 

%is(o)-^^2^<^z'^< 

X Tr[2(5^s - S^BAa^BAp - s\^pSABa], (22) 

where the scattering matrices should again be evaluated 
at the Fermi energy. The noise power of the z compo- 
nent (polarized parallel to the magnetization) of the spin 



C. Shot noise 

Shot noise of the electronic charge current is an out- 
of-cquilibrium phenomenon proportional to the current 
bias. Shot noise is due to the discreteness of the electron 
charge, and the probabilistic incidence of electrons on the 
scatterer/resistor. Let /zl — ^ir = eU with U the applied 
voltage, and take the temperature to be zero. We are here 
only concerned with the current fluctuations, although in 
this case also the average charge current is nonzero. The 
average spin current accompanying the average charge 
current does not exert a torque on a single ferromagnet, 
since the spin current is polarized along the direction of 
magnetization. From Eqs. (H^), (fTB)) . and making 

use of the zero temperature relations /^(l — /a) = and 
/ dE{fL — JrY = e|J7|, we reproduce the well-known 
charge shot noise expression^ 

^c:'!4(0) = ^|C/|[Tr(r|rTt|tT) + Tr(r|nt|tJ] (25) 

Again, the scattering matrices should be evaluated at the 
Fermi energy, and the superscript (sh) emphasizes that 

this is shot noise. 5'^'^b(0) = -S''^^2a{'^)^ where B ^ A. 
The spin current shot noise power is 

al3 CD 
X ^T:[s\caSADI3S^BDf3SBCa]- (26) 

From this we find, S^:^^j,^ = -S^:^l^ and S^^^^^^ = 

^mf^LR ~ which hold for the same reasons as for the 
thermal noise. 



D. Magnetization noise and damping 



The absorption of fluctuating transverse spin currents 
at the ferromagnet's interfaces imphes a fluctuating spin- 
transfer torque on the magnetization. The resuhing in- 
crement of the magnetization noise can be calculated us- 
ing Eq. (HI) , which by conservation of angular momentum 
is modified by the spin torque — 7oIs,abs {t) / {MsV). Here 
Is,abs = Is,L +Is,_R IS the (instantancously) absorbed spin 
current. (Recall that on both sides of the ferromagnet, 
positive current direction is defined towards the magnet.) 
Since Is.ahs is perpendicular to m, we may in general 
write ls,ahs — — m x [m x I^^abs], such that the modified 
stochastic LLG equation reads 

dm , /(!-, , rfm 

— = -7om X Hoff + h^^'it)] + a„m x — 
at at 

+ ^mx [mxl,^,j. (27) 

For the single ferromagnetic scattercr (Is,abs) = 0, but 
5Is,a.bsit) ^ 0. We can thus define h{t) = -l/{MgV)m x 
6Is (t) to be a stochastic " magnetic" field that takes into 
account the (thermal or shot) spin current noise that 
comes in addition to the intrinsic noise field h(o)(t). The 
correlators of the field 

{h,{t)h.it')) = ^^Y.S,,,AB{t~t') (28) 

^ AB 

and 

{h.{t)h,it')) = ^o^^ABit - t!) (29) 

^ AB 

for i,j = x,y; i ^ j are directly obtained from the cur- 
rent noise. h.(t) per definition has no component parallel 
to the magnetization. In the limit that the current noise 
is 'white' on the relevant energy scales (temperature, ap- 
plied voltage, and exchange splitting), we can approxi- 
mate Sij^Asit - t') w Si-j^ABi^ = 0)(5(t — t'). Using Eq. 
([22|) we then find the already advertised result 

{h'^''\t)hf\t')) = 2kBT-^5,,5{t - t'), (30) 

for the thermally (th) induced stochastic field. Here 
«' = 7o^R-e(g|-'- + g^j^)/{A-KMsV) is the spin-pumping 
enhancement of the Gilbert damping constant. This re- 
sult is in agreement with the FDT [Eq. ([3])] with a total 
Gilbert damping a = aQ-\- a' . 

Using Eq. (|26p and the unitarity of the scattering ma- 
trix we find for the stochastic field generated by the shot 
noise h'*'^) : 

(/.f \i)/.f )(i')) =^]if72%'5(^ - t')[Tr{r,r\t\t'l) 
-f T,{r[r'h^t\)] . (31) 




FIG. 2; A spin valve with two ferromagnets Fi and F2 with 
unit magnetization vectors mi and m2, here shown in the 
parallel (P) configuration mi = m2 = z. The magnetization 
of F2 is fixed. The currents in the system are evaluated by 
magnetoelectronic circuit theory on the normal side of the 
interfaces, with positive directions defined by the arrows. 

For a simple Stoner model it can be shown that for typi- 
cal experimental voltage drops in nanoscale metallic spin 
valves, h'^'''^) can dominate h*^'^^ at temperatures of the 
order of 10 K."^" In the following section we concentrate 
on room temperature, at which shot noise may be disre- 
garded. 



IV. SPIN VALVES 

We now proceed to consider the noise properties of 
spin valve nanopillars, i.e., layered structures consisting 
of two ferromagnets Fi and F2 with respective unit mag- 
netization vectors mi and m2 that are separated by a 
thin normal metal spacer N , as sketched in Fig. [51 We 
first assume that F2 is highly coercive, such that the fluc- 
tuations of it's magnetization vector are small. Such a 
'pinning' is routinely achieved in spin valves, e.g., by 'ex- 
change biasing. ' We relax this condition in Sec IIVEI 

The magnetization noise of the free layer Fi is caused 
by intrinsic processes as well as by fluctuating spin cur- 
rents in the neighbouring normal metals. The latter 
source is affected by the presence of the second ferro- 
magnet. Magnetoelectronic circuit theor y^^i'^^i^^ enables 
us to compute the current fluctuations and thus the mag- 
netizations noise of composite structures such as spin 
valves. 

Fluctuations of mi cause an easily measurable elec- 
trical noise, since the resistance of a spin valves depends 
on the relative orientation of the magnetizations (GMR) . 
Resistance noise is also interesting from a technological 
point of view, since it affects the sensitivity of spin valve 
read heads in magnetic storage devices. 

In the following, we briefly explain the spin current 
noise calculation by magnetoelectronic circuit theory. 
The stochastic field that acts on the free layer Fi and 
the related Gilbert damping are found for different mag- 
netic configurations. Using the LLG equation, we then 
calculate the fiuctuations of the magnetization vector and 
the resulting resistance noise. We finish this section by 
considering spin valves in which both ferromagnets are 
identically susceptible to fiuctuations. 
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A. Circuit theory 

Magnetoelectronic circuit theorji^i^i^S is a tool to de- 
termine transport properties of magnetoelectronic het- 
erostructures such as the spin valve shown in Fig. [21 It is 
based on the division of a given structure into resistive el- 
ements (scatterers), nodes (low resistance interconnects), 
and reservoirs (voltage sources) . The current through lo- 
cal resistors is calculated by LB scattering theory, which 
requires that nodes and reservoirs are characterized by 
(semiclassical) distribution functions. Here we take the 
ferromagnetic inserts as scatterers, the central normal 
metal layer as a node, and the outer normal metals L 
(left) and R (right) as large reservoirs. The reservoirs 
are in thermal equilibrium, and hence characterized by 
Fermi-Dirac distribution functions — f{E — fj,]^) and 
Ir = f{E — /^fl), where and fin are the respective 
chemical potentials. Depending on the relative orienta- 
tion of the magnetization vectors mi and m2, there can 
be a non-equilibrium accumulation of spins on the nor- 
mal metal node, thus characterized by a scalar (charge) 
distribution function fcN, and a vector spin distribution 
function f^^r. fcjy and fsAr form the distribution matrix 
In — ifcN + ^ ■ ^sN in 2 X 2 spin space. As before, the 
ferromagnets are thicker than Ac but thin enough such 
that spin-flip processes can be disregarded. We also as- 
sume that spin-flip in the central normal metal node is 
negligible. We are in the diffuse scattering regime, so /at 
is isotropic and constant in space. 

Referring back to Eq. (fT^ , we need now quantum sta- 
tistical averages (a\^J^E)aBnf}{E')) = 8AB&m.n5{E - 

E')f^"{E), where asnp is the annihilation operator for 
electrons moving in normal metal A {A = L, R or N) 
towards one of the ferromagnets, and is the /3q;- 
component of the 2x2 semiclassical distribution matrix 
/a in spin space. For the reservoirs (A = L ot R), we sim- 
ply have = 6f3af{E — Ha)- In contrast, in the central 
node the spin accumulation is not necessarily parallel to 
the spin quantization axis in either of the ferromagnets, 
meaning that non-diagonal (/? ^ a) terms in the distri- 
bution matrix do not vanish. The average charge current 
flowing from the right into ferromagnet -Fi can then be 
expressed by the generalized LB expressions'^^'^^ 



{Ic.ir) I dE 



glifcN + fsN • mi - /i) 



giifcN-tsN-m.^h) , (32) 



whereas the average spin current reads 

= ^ J dE |mi glifcN + isN ■ mi - /l) 
-giifcN - fsjv • mi - /l) + 2Regl\^mi x {(sn x mi) 
+2Im5j]jf,jv X mi}. (33) 

Here g" is the spin-dependent dimensionless conductance 



of Fi and gj^ is the mixing conductance of the inter- 
face between Fi and the middle normal metal. The 
average charge current and the component of the spin 
current polarized along the magnetization are conserved 
through the ferromagnet. Hence {Ic.il) — —{Ic.ir) and 
{1s,il) ■ mi = — (Is.iij) • mi. The transverse spin current 
is absorbed in the ferromagnet, leading to 

(Is,iL> = ^ J dE mi gl{fL-fcN-fsN-mi) 

-giih- fcN + tsN -mi)]. (34) 

Similar expressions hold for the currents evaluated on the 
left and right sides of F2. In order to keep the expressions 
simple we adopt from now on the parameters 5" = 52 = 
g", and gli = gl\, = g^ = g^i^ = g^ . 

Since spin-flip processes are disregarded, both charge 
and spin are conserved on the middle normal metal node: 



{Ic.ir) + {Ic.2l) = 
(Is.ir) + (I.,2l) = 



(35) 
(36) 



Eqs. (|32p - p6)) come down to four equations for the four 
unknown components of the distribution matrix /jv as a 
function of the angle 9 = cos~^(mi • m2) and the applied 
voltage U = {fiL — fJ'R)/e. Eq. ([32]) then yields {Ic il) = 
-{Ic.ir) = {Ic.il) = -{Ic,2r) = Ic = G,U, where^ 



^ e^q f , 1 — cos6' 

" 2h \ 1 — cos9 + r; + 7/cos^ 



(37) 



is the spin valve conductance with material parameters 

g = g^ +gK p = {g^ - g^)/g, and ij = 2g^^/g . 



B. Current noise 

We combine spin and charge current fluctuations, e.g., 
A/c,ifl(t) and AIs_i/j(t), respectively, on the right side of 
Fi, into a 2 X 2 matrix in spin space: 

AiMt) = iAIc.Mt) - {2e/h)d- ■ AI,4i^(^) . (38) 

Since we focus on the zero frequency noise, instantaneous 
charge and spin conservation in the central node may be 
assumed, i.e. 



A/ifl(t) + Al2L(i) =0, 



(39) 



which requires that the distribution matrix in the node 
fluctuates. The current fluctuations can then be written 

^hRi2L){t) = Shn(2L){t) + ^%^<5/^(i), (40) 



N 



where 6f]\f{t) are the fluctuations of the distribution ma- 
trix, and (5/i/j(2L) (t) are the intrinsic fluctuations [when 
'5/jv(i) = 0], coinciding with the fluctuations calculated 
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for single ferromagnets in the previous section. Expres- 
sion (|^ applies also to the current fluctuations evalu- 
ated on the left side of ferromagnet Fi and the right side 
of ferromagnet F2. In the following, we focus on thermal 
current noise, recalling from Sec. IIII Dl that for typical 
voltage drops in spin valves, shot noise is only important 
at low temperatures. 

From Eqs. (03]), ([MD and (gOD and results from 

Sec. mil we can evaluate the charge and spin current 
fluctuations in the spin valve. The correlator iS'c(O) = 
^ d{t - t'){AIc{t)AIcit')) of the charge current fluctua- 
tions is simply related to the conductance (I57|) by the 
following configuration-dependent FDT: 



(41) 



In the low-frequency regime considered here, charge cur- 
rent noise is the same anywhere in the spin valve. G„ 
can vary easily by a factor of two as a function of 9, 
which corresponds to the same variation in noise power. 
Resistance noise via magnetization fluctuations is an ad- 
ditional source of electric noise that is treated below 

The spin current correlator {AIs-^Ait)AIs-^B{t')), 
where i and j denote Cartesian components and A{B) — 
IL, IR, 2L or 2R, can be found analogously. Since spin 
current is not conserved at the ferromagnetic interfaces, 
the spin current correlator depends on the location in 
the spin valve and is not directly observable. We there- 
fore proceed to evaluate the magnetization fluctuations 
caused by the spin current noise in the next subsection. 



C. Magnetization noise and damping 

The current-induced stochastic field acting on Fi fol- 
lows from the spin current fiuctuations as explained in 
Sec. IIII Dl Here we discuss this field and, by using the 
FDT, the corresponding Gilbert damping enhancement 
in spin valves. In order to keep the algebra manageable, 
we focus on the most relevant parallel, antiparallel and 
perpendicular configurations (cos0 = O, ±1). The mix- 
ing conductances are taken to be identical for all four 
F|N-interfaces. In the semiclassical approach, intrinsic 
current fluctuations are not correlated across the node, 
implying that {SIs,,iR{t)5Is^ ,2L{t')) = 0. 



1. Parallel configuration 

For the parallel (P) magnetic configuration, mi • m.2 = 
1, the thermal spin current-induced stochastic magnetic 
field in ferromagnet Fi reads 

(/.f ) '(t'))p = 2ksT-^S,,Sit - t') (42) 

where i,j label vector components perpendicular to the 
magnetization, and 



By the FDT, asv is identical to the spin-pumping en- 
hancement of the Gilbert damping of the Fi magnetiza- 
tion. This can be checked by following the steps outlined 
for a single ferromagnet, Eqs. ([21)- ([5]). A possible ex- 
change coupling between the ferromagnets modifies the 
dynamics via Hcfr in the LLG equation, but docs not 
affect the stochastic field and Gilbert damping. 

The field correlator and damping for the parallel con- 
figuration is reduced by a factor 3/4 compared with pO|) 
for the single ferromagnet sandwiched by normal met- 
als. This result may be found also in a more direct way: 
Using Eqs. ([7]) and ([55)1 we can compute the net spin an- 
gular momentum leaving each of the ferromagnets when 
the magnetizations are slightly out of equilibrium, and by 
conservation of angular momentum infer the correspond- 
ing enhancement of the Gilbert damping constant. The 
factor 3/4 follows from the diffuse/chaotic nature of the 
node: Half of the spin current that is pumped into the 
node is reflected back and reabsorbed by Fi. 

One subtle point needs to be noted in this discussion: 
When the F-N interfaces are nearly transparent, the in- 
terfacial conductance parameters from scattering theory 
should be corrected for spurious so-called Sharvin con- 
ductances (See Sec. II. B. of Ref. i29i). In practice, this 
will correct only by a numerical prefactor close to 
one. 



2. Antiparallel configuration 
For the antiparallel (AP) configuration (mi-m2 = — 1), 

(/.f )(o/.rv'))AP = (/.r)(i)/if ^(t'))p, (44) 

i.e., the current-induced noise and damping is the same as 
in the P configuration. This result holds only when the 
imaginary part of the mixing conductance is negligibly 
small. 

3. Perpendicular configuration 

When the F2 magnetization is pinned along the x- 
direction and mi points along the z-axis 



(4"^)(t)4'h)(t'))^^2fcsr 
(/if )(t)/if)(i'))^ ^2fcsr 



a' 

yy 



5{t-t'), (45) 



6{t - 1'), (46) 



where the subscript _L emphasizes that this is valid for 
the perpendicular configuration, and, according to the 
FDT, 



^joKReg 
SttMsV 



(43) 



vv 



8ttM,V 

iTTM.V 



n 



2 - 



77(2 - P2 _|_ 2r]) 



2{l + r]){l - + rj) 



(47) 
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is the spin pumping-induced enhancement of the Gilbert 
damping. The cross correlators {hx'^\t)h^^\t')) ± = 
{hy^^\t)h^x'^\t')) ± = 0. In non-coUinear spin valves, the 
noise correlators and the Gilbert damping are therefore 
tensors. This can be accommodated by the LLG equa- 
tion for mi by a damping torque mi x Vdmi/dt, where 
the Gilbert damping tensor (in the plane perpendicular 
to the magnetization) reads: 









ao + a: 



yy 



(48) 



Note that the damping tensor must be written inside the 
cross product in the damping torque to ensure that the 
LLG equation preserves the length of the unit magneti- 
zation vector. 

In our evaluation of the Gilbert damping ((47)) , we have 
assumed that the outer left and right reservoirs have a 
fixed chemical potential which allows charge current fluc- 
tuations into the reservoirs. This is valid when the reser- 
voirs are connected to external circuit elements with suffi- 
ciently long i?C-times compared to the FMR precession 
period. In the opposite limit, when the reservoirs are 
fully decoupled from other circuit elements, charge cur- 
rent into the reservoirs must vanish at any time, and the 
chemical potentials fluctuate. This regime was consid- 
ered in Ref. 53 with the result 



yy 



SnMsV 
AttM.V 



n 



2 - 



l_p2 



(49) 



D. Resistance noise 

The fluctuations of the magnetization vector can be 
calculated by the LLG equation that incorporates the 
stochastic fields. Fluctuations in the magnetic config- 
uration affect the electrical resistance that depends on 
the dot product mi • m2. Resistance noise is an im- 
portant issue for application of spin- valve read heads i**^ 
Covington et al^ measured resistance noise in current- 
perpendicular-to-the-plane (CPP) spin valves, which are 
considered an alternative for the conventional current- 
in-the-plane spin valve read heads. We focus here on the 
zero-frequency resistance noise 



magnetization-induced noise can be the dominant contri- 
bution to the electric noise. The current densities con- 
sidered are not so high that shot noise dominates over 
Johnson-Nyquist noise, consistent with our assumption 
that shot noise may be neglected. 

In the following, we derive the resistance noise in the 
parallel, antiparallel and perpendicular configurations. 
Recall that the magnetization in ferromagnet F2 is as- 
sumed pinned. The analysis of resistance noise in the 
case of two fluctuating magnetizations is left for the next 
section. 



1. Parallel configuration 

The total stochastic field in Fi causes fluctuations 
Smi{t) = mi(t) — (mi) relative to its time- averaged 
equilibrium value (mi). For the parallel configuration 
(mi) = m2, such that the dot product of the magne- 
tizations is cos 6* = mi • m.2 = 1 — 6nil/2, with 9 the 
angle between the magnetization directions. For small 
fiuctuations we can expand the resistance to first order 
in ^mi 



such that the resistance noise correlator becomes 



(51) 



{AR{t)AR{t'))p = (R{t)R{t'))p - {R{t))p{R{t'))p 

-{Snilit))p{Smlit'))p] , (52) 

where the brackets denote statistical averaging around 
the parallel configuration. Assuming that the stochastic 
fields are Gaussian distributed, so are the fluctuations of 
the magnetization vectors, since the magnetization is a 
linear function of the stochastic fields. We may then em- 
ploy Wick's theorem,**^ according to which fourth order 
moments of the fluctuations can be expressed in terms of 
the sum of products of second order moments. We then 
arrive at 



{AR{t)AR{t'))p 



1 r oRji 

2 \9cos 



2 



(53) 



Sr{uj^O)^ d{t-t'){AR{t)AR{t')), (50) 



where AR{t) is the time dependent deviation of the re- 
sistance from the time-averaged value. 

Resistance noise can be measured e.g. as voltage noise 
for constant current bias or as current noise for a constant 
voltage bias. The resistance noise comes on top of the 
Johnson-Nyquist noise discussed in Sec IIVBI and in Ref. 
l43i . We find that at relatively high current densities, the 



where i and j denote Cartesian components. 
Eq. JST]) we find 



dRjl) 
9 cos 6' 



hP' 



From 



(54) 



Since the magnetization fluctuations are small, we may 
disregard their longitudinal component, whereas the cor- 
relator of the transverse fluctuations can be computed by 
the LLG equation. 
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We use the coordinate system in Fig. [2]with interfaces 
in the xz-plane. The LLG equation reads 



dm.1 
dt 



7omi X [Heff + h(i)] 



dm.1 
"dT' 



(55) 



where the total stochastic field h(i) = h^^) {t)+h^*-^'> {t) in- 
cludes both the intrinsic field h^^^^ {t) (see section lllB and 
the current induced field h'^^^^t) from the previous sec- 
tion, ao and asv are the corresponding Gilbert damping 
parameters. The effective field Hos = Ho+Ha+Hd+Hg 
contains the external field Hq, the in-plane anisotropy 
field Hq, the out-of-plane demagnetizing field H^, and 
the sum of dipolar and exchange fields Hg. The external 
and anisotropy fields are both taken along the z-axis. We 
parametrize these fields by loq and uja as 7H0 — ujqz and 
jUa — Wa(mi • z)z. The demagnetizing field is directed 
normal to the plane, i.e. along the y-axis, such that 
jUd = — ajd(mi • y)y thereby introducing the parameter 
ujd. The dipolar and exchange couplings are described 
in terms of a Heisenberg coupling — Jmi • m2, which fa- 
vors a parallel magnetic configuration for J > and an 
antiparallel one for J < 0. This translates into the field 
7He = uJeUi2, whcrc uje — 'yJ/M^d. 

In the P configuration (mi) is aligned with the pinned 
ni2 in the +z direction, which can always be enforced by 
a sufficiently strong external field. Linearizing the LLG 
equation in the amplitude of the transverse fluctuations 
dui{t) w d'mx{t)x + 5my{t)y, we find the magnetization 
noise correlator 



{5mi{t)5mj{t')) p 



dwe-''^(*-*''[/,j, (56) 



ttMsV 

by using the correlators of the stochastic fields. Here 



Uxx 

u. 



-iLu{2uJt + LUd) 



[uj^ - iutiiot + iUdW + io^a^{2iOt + WdY 

(..2+^2) 



, (57) 
, (58) 
, (59) 



Uyx - -U^y, (60) 

with a = ao + o-sv and LOt = wo + -I- Wg. The above 
expressions hold for small damping, i.e., ao,a^^ ^ 1. 
The zero- frequency resistance noise 5'p(0) = J d{t — 
t'){AR{t)ARlt'))p is obtained by inserting Eq. ^ into 
Eq. (ET 



Sp{0) =- 

TT 



1 /hp^ 



duj{a 



2 

XX 



u. 



yy 



2C/, 



xy) 



(61) 



demagnetizing field, which serves to stabilize the magne- 
tization in the plane of the film, is important to get the 
right magnitude of the noise, we can gain physical un- 
derstanding by disregarding it. Setting ujd — Q, we find 



Sp{0) 



(62) 

<a 



Obviously, the resistance noise strongly depends on the 
parameter Wf. The external and anisotropy fields sta- 
bilize the magnetization, hence lowering the noise. The 
dipolar and exchange field either stabilizes or destabilizes 
the magnetization, depending on the sign of the coupling 
constant J. We observe that the Gilbert damping also 
strongly affects the resistance noise. The resistance noise 
decreases with increasing damping, because the suppres- 
sion of the magnetic susceptibility by a large alpha turns 
out to be more important than the FDT-motivated in- 
crease of the stochastic field noise. Since can be of 
the same order as ao^^, the importance of spin current 
noise and spin pumping is evident. 

When a constant voltage bias is applied, the resis- 
tance noise causes current noise. At sufficiently small 
bias, the Johnson-Nyquist current noise (Sec. lIVBp al- 
ways wins. However, at relatively high current densi- 
ties, the effects of the resistance noise are very signifi- 
cant. That noise may be important for the next gen- 
eration magnetoresistive spin valve read heads 4^ For a 
quantitative comparison, which depends on many mate- 
rial parameters, it is important to use Eq. ()61|) and not 
Eq. (|62p . since the demagnetizing field has a large effect 
on the magnitude of the magnetization-induced noise. 
The magnetization-induced noise is most prominent for 
small structures, since the ratio of Johnson-Nyquist noise 
to magnetization-induced noise scales with the volume of 
the ferromagnet. 



2. Antiparallel configuration 

When J < 0, the dipolar and exchange coupling fa- 
vors an AP configuration ((mi) = — m2) at zero exter- 
nal magnetic field. Following the recipe of the previous 
subsection, we find a resistance noise 



{AR{t)AR{t'))AP ^^^^ 



2 V 9cos( 



x^((5mi,,(i)<5mi,,(i')>AP, (63) 



where the sensitivity of the resistance to the fluctuations 
is 



dRj-l) 

5 cos 6* ^ " 62.9(1 - P2)2- 



(64) 



To gain insight into this rather complicated expression, it 
is convenient to make some simplifications. Although the 



Using the magnetization noise correlators from the lin- 
earized Eq. (j55p . the zero frequency resistance noise be- 
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comes 



Sap{0)^ / d{t-t'){AR{t)AR{t')) 



AP 



TT Ve2c/(1-P2)2 J y ]^,J^Y 

X / dc.(14l + - 2V,%), (65) 



where Vij = Uij{uJt ~^ ^s) with uj^ — uja — tOe (recall 
that We < 0). Again disregarding the demagnetizing field 
strongly simplifies the expression: 



Sap{0) - 



MsV 



e^g{l-P^y 



(66) 



As expected, the resistance noise decreases with increas- 
ing uJs- The anisotropy, dipolar and exchange fields sta- 
bilizes the magnetization, playing a role similar to that 
of the external field in the P configuration. The Gilbert 
damping enters in the same way as for the P configura- 
tion. 

Except for the prefactor that reflects the sensitivity of 
the resistance to the magnetization fluctuations, S'p(O) 
and SApiO) are very similar. For the special case ujt — i^s, 



Sap 



{l~P^ 



(67) 



For, e.g., P = 0.7 and r] — 1, this becomes Sp/Sap ~ 6% 
showing that the difference in noise level between the P 
and AP configurations can be substantial. 

This asymmetry in the noise level between the P and 
AP configurations is consistent with the the experimen- 
tal results of Covington et al. on nearly cylindrical mul- 
tilayer pillars42. In these experiments the magnetizations 



were aligned parallel when the external magnetic field 
reached about 1500 Oe. Although we treat spin valves 
with two ferromagnetic films and Covington et al. dealt 
with multilayers of 4-15 magnetic films, it is likely that 
the difference between the noise properties of bilayers and 
multilayers is small, as the only local structural difference 
is the number of neighboring ferromagnets. This asser- 
tion is supported by the experiments by Covington et 
al. that did not reveal strong differences for nanopillars 
ranging from 4-15 layers. 



3. Perpendicular configuration 

We now investigate the perpendicular state (mi) •m2 = 
0, assuming that m2 now has been pinned in the x- 
direction, whereas mi is on average parallel to the z axis, 
as before. In the following we assume that the interlayer 
exchange and dipolar coupling are negligibly small, since 
otherwise the algebra and expressions become awkward. 

Expanding the resistance to first order in the fluctua- 
tions (5mi, we find in this case 

{AR{t)AR{t'))^ = fl^) {Sm,,it)Sm,^,{t')). (68) 
\ a cos J 

The magnetization fluctuations affect the resistance noise 
in the perpendicular configuration to second order, unlike 
for the P and AP configurations, in which the leading 
term was of fourth order. The sensitivity of the resistance 
for this configuration is according to Eq. |37|) 



dRjO) 
dcosO 



AhP'^-q 



(69) 



Linearizing Eq. (|55p and using the correlators Eqs. ([l5 
and (l46t for the stochastic field we find 



ttMsV 



dujt 



jit-t')_ 



a;2(ao -I- a'yy) + {ujp + ujd)'^{ao 



[uj^ - ujp{ujp + ujdW + w2[u;p(2ao + <^ -I- a' ) + u}d{ao + <^)]2 



(70) 



where tOp = luq+i-^c 
resistance noise 



We then arrive at the zero- frequency 



ihP'^r] 



MsV ye^gil + Tj - P^) 



(71) 



quite different from that in the coUincar configurations. 
In particular, the damping appears here in the numerator 
and there is no dependence on the demagnetizing field. 
Notice that since S± is quadratic in magnetic fluctuations 
[see Eq. (j68p ]. it becomes linear in temperature, unlike Sp 
and Sap- 



E. Two identical ferromagnets 

We now investigate spin valves in which the ferro- 
magnets are identical and hence equally susceptible to 
fluctuations,''^ focusing now only on the P and AP con- 
figurations. The fluctuations of Fi are dmi{t) — mi{t) — 
(mi) and those of F2 are 6m2{t) = m2{t) — {^2)- As 
before, we choose the z-axis so that the time-averaged 
equilibrium values are (mi) — (m2) = z for the par- 
allel configuration, and (mi) = — (m2) — z for the 
antiparallel. The dot product of the magnetizations is 
mi ■ m2 = ±1 =F ('5m^)^/2, where the upper (lower) sign 
holds for the P (AP) orientation and Sm^ = dmi ^6m2. 
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For small fluctuations, we can expand the resistance to The resistance noise is then 
first order in (^m^)^, finding 

ii(m,.m.)«i?(±)^(.m^f^. (72) 



{AR{t)AR{t'))p^AP = {R{t)R{t'))p,AP - {R{t))p/Ap{R{t'))p/AP 

= i (^S^)' W^'^?{5m^f)p/AP {{5vn^f)p,Ap{{5rn^f)p/Ap] , 



(73) 



r 



which by employing Wick's theorem becomes 



{AR{t)AR{t'))p/AP 



_ 1 ( dR{±l) Y 



P/AP- 



(74) 



Letting the subscripts k and I refer to ferromagnet 1 
or 2, the LLG equation in this case reads 



dt 



'7omfe X [Hcff + hfe(t)] 
dm.k , 



+ (ao + asv)m.k x 



dm; 



dt 



-m; x 



dt 



(75) 



where the effective field Hcff is now taken to be equal 
for both ferromagnets. Due to current conservation, the 
ferromagnets respective current-induced stochastic fields 
are not independent of each other. With the spin current 
noise calculated in Sec. IIVBI and following the recipe in 
Sec. HiTdI we find 



for the P configuration, and 



2kBT "7/i %(5(t - t') (76) 



jAP 



2ksT^^6,,S{t~t'). (77) 



for the AP configuration (as before i,j label components 
perpendicular to the magnetization direction) . a^v is de- 
fined in Eq. dig). Naturally, the bulk fields h^"^ and h^"^ 
are uncorrelated. The last term in the LLG Eq. ([75]) 
represent the dynamic spin-exchange couplingj^Si^ It is 
the spin current pumped from ferromagnet I (see Sec. lIIip 



that is transmitted to and subsequently absorbed by fer- 
romagnet k. Since the normal metal node is chaotic, 
this amounts to one third of the net total spin current 
pumped out of ferromagnet I. This dynamic coupling was 
not present in spin valves in which one magnetization is 
not moving at all. 

By linearizing Eq. ([75)1 in 6mk{t) we can evaluate the 
desired magnetization noise correlators that are to be 
inserted in Eq. (|74p . The zero-frequency resistance noise 
for the P and AP configurations then respectively reads 



SpiO) 



hP^Y f2-fokBT 



X I duj{Zl^ + Zly-2Zly). 



and 



Sap{Q) = 



1 



j duj{Xl+X, 



' ' 2^okBT ^ ' 



Here 
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\up- - Wi(tJi ■ 




u}'^al{2uji - 






4u!at{2uji 


+ Wd) 




\up- - LUt{uJt - 


+ OJdW + 


w^Qfj (2lJj - 






af(w2 + 






[UJ^ - LUi{lOi - 


+ OJdW + 


LL!'^al{2ix}i - 


Vujdf 



and 



(78) 



(79) 

(80) 
(81) 
(82) 



+ (cj^ + UJd){2uJe - iOc)Y + l-0'^(2ijJxOLs - 2uJ(,OL - UJdOLtY ' 
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FIG. 3: The resistance noise in the P configuration as a 
fiinction of the externally applied magnetic field, given in 
units of {IQ-"^ /■K){hP^ /e^gr)f{2'yokBT/MsVf. The parame- 
ters used are qq = Q-sv = 0.01, uoc/^o = i^dHo ~ 100 Oe and 
J = -0.10 erg/cm^ 

For convenience, we defined as — ao + 2as„/3, at = 
ao + Aasv/-^, a = ao + asv (note the difference between 
a, as and at), and uji = ujq + LUa + '2ujx. The above 
expressions hold for small damping, i.e., aQ,a'^y ^ 1. 

Compared to the results in the previous section, we see 
that Eq. ^ is similar to Eq. jH]), whereas Eq. ([79]) 
differs considerably from Eq. ([55)1 . This is due to the 
static dipolar and exchange couplings, and the dynamic 
spin-exchange coupling, whose effects on the noise are 
modified by the presence of the second fluctuating fer- 
romagnet. In particular, the latter coupling causes the 
Gilbert damping constant to enter Eqs. ([78]) and ([79]) dif- 
ferently. Eq. ([75)1 decreases with the external field and 
Eq. (|75)) decreases with the dipolar and exchange cou- 
pling, as expected, and as shown in Figs. [3] and ID The 
noise level is in general higher when both fcrromagncts 
fluctuate, than when only one does. 

The resistance noise is governed by a number of ma- 
terial parameters. Depending on these parameters, the 
noise level in the P configuration can differ substantially 
from that in the AP configuration. Note that Eq. ([75)1 
reduces to that of Ref. 36 when the demagnetizing field 
is disregarded, i.e., when ujd 0, whereas Eq. ([78]) 
does when Wo — > and ujd — > —(^a, since the exter- 
nal field in our earlier work was perpendicular to the 
anisotropy field. The considerable difference between 
Sp{0) can S'ap(O) in typical experimental spin- valve se- 
tups can partly be explained by the dynamic exchange 



coupling .^^ However, also the sensitivity of the resistance 
to magnetic configuration changes can be important, as 
shown in the previous section. The demagnetizing field 
also significantly affects the numerical result for the noise 
level since it stabilizes the magnetization, in both the P 
and AP configurations. 
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FIG. 4: The resistance noise in the AP configuration as 
a function of the dipolar and exchange coupling between 
the ferromagnets, given in units of (10"V7r)[/iP^»?/e^g(l - 
P^ff{2-fokBT/MsVf. The parameters used are ao = a^v = 
0.01, ajc/70 = t^d/7o = 100 Oe, and Ho = 0. 



V. CONCLUSIONS 

Using scattering theory and magnetoelectronic circuit 
theory, we demonstrate the effect of spin current fluctua- 
tions on the magnetization in ferromagnetic multilayers. 
Via a fluctuating spin-transfer torque, the current noise 
causes signiflcantly enhanced magnetization noise, which 
in spin valves is a function of the magnetic conflguration. 
The noise is related to the magnetization damping by the 
FDT, and can be experimentally detected as resistance 
noise. The contribution from spin current noise to resis- 
tance noise is considerable, and may be an issue for the 
next generation magnetoresistive spin valve read heads. 
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